SIGNATURE SIMULATOR 
Background 

This application claims priority under 35 U.S.C. § 
119(e) (1) to co-pending U.S. Provisional Patent Application No. 
5 60/428,205, filed November 20, 2002, and entitled "Methodology 
to Synthesize Agent, Simulant, Interferent or Toxic Industrial 
Chemical Signatures for an Interferometer-Based Spectral 
Sensor" , wherein such document is incorporated herein by 
reference . 

10 The U.S. Government may have certain rights in the present 

invention. 

The present invention relates to simulation and 
particularly to simulation of scenarios for testing detectors 
and related systems. More particularly, it relates to scenarios 

15 that include simulated interferometer-based signatures of agent, 
stimulant, interferent and toxic chemicals. 

Testing of detectors or systems in real situations and 
scenarios involving field measurements may often be very 
expensive and quite hazardous. In lieu of the field 

20 measurements, traditional simulations have been used. However, 
they have been often inadequate for testing detectors, sensors 
and related systems. Typically, traditional simulations have 
been predicated on the principles of physics and do not 
necessarily represent accurately the phenomenological behavior 
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needed for the testing the detectors, sensors, systems, and the 
like . 

Summary 

The present invention is an approach for constructing, 
calculating and/or synthesizing agent, simulant, interferent, 
toxic or non-toxic chemical signatures to provide possible 
scenarios for testing an interferometer-based spectral sensor, 
detector, system, or the like. The simulations may address the 
artifacts and parameters specific to a distinct sensor, detector 
or system. The present signature synthesis methodology may be 
different from that in the related art because of its relatively 
precise representation of sensor artifacts, sensor sampling 
error and sensor signal-to-noise ratio. 

Brief Description of the Drawing 
Figure 1 shows a block diagram of an algorithm design 
me t hodo 1 ogy ; 

Figure 2 shows the geometry of the two- slab cloud model for 
the agent/ simulant/interf erent simulation; 

Figure 3 is a block diagram of a baseline simulation; 

Figure 4 reveals the computational stages of the baseline 
simulation; 
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4. 

Figure 5 is a top-level diagram of a detector chemical 
agent simulation; 

Figure 6 is a diagram of the preparation stage for the 
simulation; 

5 Figure 7 is a diagram of the calibration stage of the 

simulation; 

Figure 8 is a diagram of the restoration stage of the 
simulation; 

Figure 9 is a diagram of the major components of the 
10 simulation software; 

Figure 10 shows a profile of the values from a start scan 
to an end scan in a simulation run; 

Figure 11 high level operational flow diagram of the 
simulation software; and 
15 Figure 12 shows the removal of the original apodization 

effect by dividing the truncated interf erogram with a window. 

Description 

The invention may be an interferometer-based chemical agent 
20 simulator which may synthesize agent, simulant, interf erent 
and/or toxic or non- toxic chemical signatures in background 
interf erograms . The signatures may be collected by an arbitrary 

3 
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interferometer sensor. Further, there may be simulated sensor 
outputs that have similar characteristics of a distinct, 
possibly non-existent sensor. The simulator may produce 
simulations to be integrated into an automated experiment 
5 configuration and/or performance evaluation environment to test 
sensors, detectors, and related systems. One detector, sensor 
or system that may be tested is a Joint Services lightweight 
stand-off chemical agent detector (JSLSCAD) . The JSLSCAD may be 
regarded as a state-of-the-art interferometer-based chemical 

10 agent detector having baseline chemical agent detection 

algorithms. Although the present simulator may be used for 
testing many kinds of sensors, detectors and systems, the 
JSLSCAD may be considered as an example detector that may be 
provided simulated chemical agent signatures. 

15 The simulator may be designed with algorithms to calculate 

and/or synthesize artificial interferometer-based signatures for 
a simulation of a scenario to be sensed. Figure 1 shows a block 
diagram of an algorithm design methodology. Field data 11 and 
atmospheric conditions 12 may be input into a ground truth tool 

20 13 . Ground truth tool 13 may allow the user to enter or 

establish the ground truth, e.g., types of background, and tool 
13 may compute some basic features of the input, e.g., peak 

4 



H0003 7 98 (1100 .1204101) 

amplitude, locations and kurtosis. (Kurtosis may be a quantity 
indicative of a general form of a statistical frequency curve 
near the mean of a distribution.) Selected data with ground 
truth from tool 13 may go to a simulation block 14. Also input 
5 to simulation block 14 may be agent absorption coefficients 15 
and system performance requirements 16. Simulation 14 may add 
chemical agent and/or common battle field signatures to 
background data which cannot be obtained in field collection. 
Simulation 14 may have training data 17 output and a sequestered 

10 test data 18 output. Sequestered test data 18 may be input to 
detection performance test and evaluation 19. Training data 17 
may be input to algorithm 20. An output of algorithm 20 may go 
to detection performance test and evaluation 19. Pd 
(probability of detection) and Pfa (probability of false alarms) 

15 signals may be an output of Detection performance test and 

evaluation 19. This output may be fed to algorithm development 
20. Detection performance test and evaluation 19 may test the 
sensor or detector outputs (both field and simulated data) with 
a suite of algorithms. It may also evaluate the performances 

20 via a comparison with ground truth information 13 . Algorithm 
development 20 may design, code and train detection algorithm 
candidates based on the sensor, agent/interf erents and 

5 
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background phenomenologies . Features of the algorithm design 
methodology may include model based on first principle 
phenomenology and a unified simulation, development /train, test 
and performance evaluation package. 

A baseline detector or sensor chemical agent simulation 25 
may be designed for calculation of synthetic chemical agent 
signatures utilizing background scene measurements obtained with 
an arbitrary sensor 21 of Figure 2. A measurement may be 
transformed into a synthetic representation of the sensor data 
that would be obtained with a possibly distinct sensor, 
utilizing chemical agent absorptivity coefficients, 
specifications for the sensor's sample comb, the scene 
configuration and the signal-to-noise ratio of sensor 21. The 
baseline simulation configuration for signature synthesis is 
shown in Figure 2 with radiance (N C i) and transmission (x C i) of 
the front cloud slab CI, radiance (N C2 ) and transmission (x C2 ) of 
the back cloud slab C2 , the background of the scene with 
radiance (N B ) and foreground (N F , t f ) . 

A block diagram of simulation 25 is shown in Figure 3. The 
major components of simulation 2 5 may include sensor response 
removal 24, spectral characteristics addition 27, atmospheric 
attenuation 28, MODerate-resolution atmospheric TRANmittance and 
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w 

radiance model (MODTRAN) 23, and sensor response addition 29. 
" MODTRAN is a FORTRAN- like computer code designed to determine or 
model atmospheric transmission and radiance. MODTRAN was 
developed principally by the U.S. Air Force Research 
5 Laboratories. Simulation 25 may include sensor response removal 
having inputs from field data 11 and sensor response 26, and an 
output connected to an input of a spectral characteristics 
addition 27. Agent absorption coefficients 15 may be connected to 
another input of special characteristics addition 27. An output 

10 of special characteristics addition 27 may be connected to an 
input of atmospheric attenuation 28. MODTRAN 23 may have an 
input from atmospheric conditions 12 and an output pertaining to 
radiance and transmission of the atmosphere connected to another 
input of atmospheric attenuation 28. An output from atmospheric 

15 attenuation may be connected to an input of sensor response 
addition 29. Sensor response 26 may be connected to another 
input of sensor response addition 29. An output of sensor 
response addition 2 9 may be simulated data. 

Sensor response removal 24 may remove the sensor induced 

20 affects in the input signal and the output is sensor- independent 
signature of the background. Spectral characteristic addition 
27 may have the signatures of the agent and /or battlefield 

7 
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interferents added to the background only signature. The 
signatures may be based on the concentration, temperature 
difference, and field of view (FOV) coverage. Atmospheric 
attenuation 28 may have the signature properly adjusted due to 
the atmospheric transmission, based on the range of the agent 
cloud and atmospheric conditions. Sensor response addition may 
have the affects of the sensor response introduced into the 
signature, creating the simulated data that the sensor should 
output . 

Feature's of simulation 25 may include a first principle 
phenomenological model (spectral absorption/radiation) and 
Government off the shelf (GOTS) proven atmospheric models may be 
used relative to MODTRAN 23. Simulation 25 may be sensor 
independent (i.e., may simulate signature of arbitrary Fourier 
transform infrared (FTIR) spectrometers, including JSLSCAD) 
Simulation 25 may have two agent /interferent clouds that allow 
multiple agents and battlefield entries, which can be either 
mixed or separated. 

The computational stages of the baseline simulation 25 for 
the detector are shown in figure 4 . Parameters for each 
spectrum 46 may be computed. From ground truth of tool 13 and 
from case definition, specific parameters may be determined for 

8 
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each subfile. These parameters may include FOV coverage 
parameters, MODTRAN 23 inputs, cloud slab (CI and C2) 
temperatures, components and their CLs . A calibrated background 
spectrum 51 (N B ) may be obtained. Relative to an atmospheric 
model 52 from conditions 12, a MODTRAN 23 input file may be 
constructed. MODTRAN 2 3 may be run and the radiance N F and 
transmission x F arrays may be recovered. N F and x F may be 
convolved with a sine function representing the Fourier 
transform of the temporal averaging window corresponding to the 
sample spacing of the interferometer to adjust resolution. 
Also, N F and x F may be interpolated to sample a comb of N B . 

A cloud model 53 involving N C i, t C i and N C 2 and t C 2 may be 
determined. For each cloud slab (CI and C2) , an argument (ARG) 
may be constructed for the components of the respective slab. 

ARGci = (XiCLi + a 2 CL 2 + . . . + a n CL n 
ARG C2 = OCiCLi + a 2 CL 2 + . . - + a^CL™ 

where oti = absorptivity coefficient for component i as a function 

of v. Transmission terms x C i = exp ( -ARG C i) and x C2 = exp ( -ARG C2 ) 

may be computed. Also computed may be radiance terms N C i = 

N BB / ci (1-Tci) and N C2 = N BB , C 2 (l-t C 2) where N BB/C = Planck's function 

for each slab at the temperature of the slab as a function of v. 
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A simulated spectrum 47 (N S imul) may be built as indicated by 
the following equation. 

Nsimul = Si* (N B -N F ) *Xci*Tc2 + Si*Nc2*T F *Tci + S 2 * (N B - 

N P ) *Tc2 + S 2 *N C2 *t f +Si*N C i*Tf+N p , 

where Si = fraction of FOV coverage by cloud 1 and S 2 = l-Si. Then 

a signal-to-noise compensation 48 may be incorporated. 

An alternate realization of the phenomenological synthesis 
calculation may be the following. If the sensor-based 
measurement is defined to be the background of the scene, then a 
methodology may exist to estimate the signature that would be 
observed if the sensor was actually further away from the scene 
along the same view direction. The reference scene does not 
have to be compensated for the atmosphere between the sensor and 
the agent cloud because the agent cloud may be envisioned to be 
located at the sensor location. The "virtual" sensor may be 
further away from the scene along the view direction of the 
sensor for the original sensor data. Then the equation 
designated "Build simulated spectrum" in Figure 4 may become 

Nsimul = Si*N b *Tci*Tc 2 *Xf+Si*Nc 2 *t f *Tci + s 2 *Nb*x C2 *t f + 

S 2 *N C2 *T f +Si*Nci*T f +N f , 

where 

Si = fraction of FOV coverage by cloud 1 and 

10 
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y 

s 2 = 1-Si 

Figure 5 is a top-level diagram of the detector (e.g., 
JSLSCAD) chemical agent simulation. It shows inputs 31 (to 
simulation 25), preparation 32, calibration 33, simulation 25 
5 and restoration 35. Inputs 31 to simulation 25 may include a 
background interf erogram and ground truth, ambient and liquid 
nitrogen blackbody interf erograms and ground truth, absorptivity 
data and a case definition of the simulation. Preparation 32 
may include preprocessing 36, data conversion 37, normalization 

10 38 and fast Fourier transform (FFT) 39. Calibration 33 may 
include compute ambient blackbody spectrum 41, compute 
theoretical ambient blackbody spectrum 42, compute theoretical 
LN 2 spectrum 43, compute response correction array 44 and compute 
calibrated spectrum 45. Simulation 25 may include compute 

15 parameters for each spectrum 46, build simulated spectrum 47 and 
signal-to-noise compensation 48. Build simulated spectrum 47 
may have inputs background spectrum 51 (calibrated N B ) , 
atmospheric model 52 (N F , t f ) , and cloud model 53 (N C i, t C i, N C 2 / 
t C 2) . Restoration 35 may include uncalibrate 56, inverse fast 

20 Fourier transform (IFFT) 57 and output to database 58. Figure 4 
is a simulation stage 25 block diagram. 
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Figure 6 is a preparation stage 32 block diagram. 
Preprocessing 36 may select a set of subfiles upon which to 
simulate agent/interf erent signatures (allowing for co-adding, 
if desired) and perform co-adding, if desired, to yield a 
5 contiguous set of subfiles to utilize. Data conversion 37 may 
transform measured interf erogram signatures to sample spacing of 
a new interf erogram sensor. Normalization 38 may normalize for 
sensor gain and, if co-adding is performed, the data should be 
normalized (i.e., data points/# of co-adds per resultant 

10 subfile) . FFT 3 9 (I B ) may compensate sensor measurement for 

sensor self radiance (pointwise subtraction of a liquid nitrogen 
calibration interf erogram from a measured interf erogram) , 
triangular apodization (i.e. a tapering function), zero pad to 
4096, FFT and phase correct and preserve reference phase for 

15 restoration process. 

Figure 7 is a calibration stage 33 block diagram. Compute 
ambient blackbody spectrum 41 (I A ) may include a pointwise 
average of subfiles 1-128 from calibration file and a 
compensation for sensor self radiance, triangular apodization, 

20 FFT and phase correct. Compute theoretical ambient blackbody 

spectrum 42 (N bb ,a) may be at T = Tambient + 273.15 degrees Kelvin, 
calculated from a Planck function for an I A sample comb. Compute 

12 
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theoretical LN 2 spectrum 43 (N bb ,l) may be at T = 77 degrees 
Kelvin, calculated from a Planck function for an I A sample comb. 
Compute response correction array 44 (R) where R may be 
calculated as I A / (N BB#A -N BB/L ) . Compute calibrated spectrum 45 (N B ) 
where N B may be calculated as I B /R- 

Figure 8 is a restoration stage 35 block diagram. 
Uncalibrate 56 (I S imul) where I S imul may be calculated as N S imul*R. 
IFFT 57 may subtract the original background spectrum, use a 
stored reference phase from FFT to reapply phase, IFFT, remove 
the apodization originally applied, reintroduce sensor self- 
radiation and unnormalize for gain. Output to database 58 may 
output simulation results to a database properly annotated with 
simulation ground truth parameters. 

There may be simulation software that adds the chemical 
agents and/or common battlefield inferents into background 
interf erograms . Simulation software may synthesize agent, 
simulant and/or inferent signatures in background interf erograms 
which may be collected by a sensor. The outputs may be sensor 
specific and expected to have similar characteristics of the 
anticipated sensor. The simulation may address many artifacts 
and parameters specific to the sensor. A simulation that may 
apply to a different sensor would require modification to the 

13 
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simulation software. The terms "scan" and "interf erogram" may 
be used interchangeably. A scan or interf erogram may mean a 
2048 point data sequence that is or simulates the sensor output. 
Each data point of the sequence may be represented by a 16 bit 
word. 

The assumptions and requirements needed to complete the 
simulation software may include a local, executable MODTRAN 

(version 3.7) which can be invoked by the simulation software, 
an executable MATLAB® (version 5.2) and its signal processing 
toolbox which can be invoked by the simulation software. MATLAB® 

(Matlab) may be used as a numerical computing tool. The 
background interf erograms used for simulation may be already 
imported in chemical agent detection software environment 

(CADSE) database. Also, these interf erograms may be already 
converted to have the same attributes as the JSLSCAD, e.g., 4 
wavenumber resolution and 2048 points. The ancillary 
information associated with the interf erograms (e.g., ambient 
blackbody and liquid nitrogen references) may be also correctly 
stored in the CADSE database. The absorption coefficients of 
the agent and interf erents may be stored in a file whose format, 
resolution, and number of points are specified. 



14 
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The simulation software may be partitioned into three major 
• components: CADSE simulation 61, Matlab simulation 62 and 

MODTRAN 23. A user interface 63 may be part of CADSE simulation 
61. Figure 9 shows the connections between the components. 
5 User interface 63 may have a two-way connection with CADSE 

simulation 61. CADSE simulation 61 may have input connections 
from atmospheric transmittance and radiance 65 , background 
measurement experiment or environment (BME) interf erogram 66 and 
ancillary information 67. CADSE simulation 61 may also have 

10 two-connections with MODTRAN 23 and Matlab simulation 62. CADSE 
simulation 61 may have output connections to simulated 
interf erogram 64 and cards 68. MODTRAN 23 may have an input 
connection from cards 68 and an output connection to atmospheric 
transmittance and radiance 65 . Atmospheric transmittance and 

15 radiance 65 may be the atmospheric information from MODTRAN 23 
to CADSE simulation 61. 

CADSE simulation 61 may be the coordinator that obtains the 
necessary information from the user via user interface 63, gets 
BME background interf erogram 66 and associated or ancillary 

20 information 67 from the CADSE database, creates the protocols to 
MODTRAN 23 and Matlab 62, and saves simulated interf erograms 64 
and their ground truth information into the CADSE database. 
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MODTRAN 2 3 may be the software for standard atmospheric 
» transmission and radiance 65 model. MODTRAN 23 may be used to 
produce foreground atmospheric transmission and radiance. 
Matlab simulation 62 may generate the simulated 
5 interf erograms 64 based on a two-slab cloud model. All of the 
mathematical computations may be performed in Matlab simulation 
module 62 . However, the simulation software may run using 
Matlab or C. The simulation software may require numerous 
inputs, which can be grouped into three categories: MODTRAN 23 
10 related, background interf erogram related, and cloud slab 

related. All the input values for a simulation run may be saved 
and recalled for future reference. 

The following items pertain to MODTRAN 23 related 
parameters. Sets of eight cards 68 (i.e. input files in card 
15 format) may be utilized for providing required inputs to the 
MODTRAN 23 model. The values of most parameters may be fixed 
for this application. Values of some parameters may vary 
depending on the scenario, which the user may input via 
interface 63 . Alternative values may indicate options that are 
20 related to different background and simulation scenarios. CADSE 
simulation 61 may develop the sets of cards 68 based on the user 
inputs and the above-noted fixed parameter values. These sets 

16 



H0003 7 98 (1100 . 12 04101) 

of cards 68 may be saved as a file that is input to the MODTRAN 
23 executable. 

The background interf erogram related parameters may be a 
category of inputs that relate to the background interf erograms 
66 in which agent and/or interf erent signatures are to be added. 
Some inputs may specify the background interf erograms 66, and 
some inputs may relate to the background interf erograms 66, such 
as the calibration references. The majority of these inputs may 
be retrieved from the CADSE database once they are specified. 
The list of inputs, their format and descriptions are tabulated 
in Table 1. 



Table 1: Input parameters related to the input interf erograms 66 



Name 


Format 


Value/ 
source 


Description 


Scan set 


CADSE 
UI 




The scan set specifies the set of interferograms in 
which agent/interferents will be simulated upon. 


Ambient 

Blackbody 

Reference 




Calibrati 
on file 
referred 
in global 
file 

header 


Each file that a selected scan originated from is 
associated with a calibration file. The calibration file 
consists of three sets of references: an ambient 
blackbody set, a polystyrene reference, and a LN 2 
reference set. The ambient blackbody set that 
corresponds to the selected scan is needed. 


LN 2 Reference 




Ditto 


This is the LN 2 reference set as stated above 


Slant angle 




Subfile 
header 


The slant angle of the scan. Noted that the angle 
required for MODTRAN is Zenith angle, not the 
elevation angle as in the subfile header. 


Gain 




Global 
header 


The gain of the sensor. 


Ambient 
temperature 




Global 
header 


The ambient temperature when the scan was collected. 


Barometric 




Global 


The barometric pressure when the scan was collected. 
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Pressure 




header 




Relative 
humidity 




Global 
header 


The percentage of relative humidity when the scan was 
collected. 


numCoadd 




User 


Number of interferograms coadded 



A byproduct of the input information may be the number of scans 
(subfiles) which the agent and interferents will be simulated 
upon. This number may be referred as numScan throughout this 
document . 

The following pertains to cloud slab related parameters. 
Two cloud slabs may be simulated in the model. Further, the 
model may allow multiple components in each cloud slab. Here, 
the maximum number of components may be four, even though the 
number may other than four. The model also may assume that the 
back slab has 10 0 percent FOV coverage, and that the front slab 
may have 0 to 100 percent FOV coverage. To describe the clouds 
of agent/interf erents , following inputs listed in Table 2 may be 
considered to be required. 



Table 2 : Input parameters related to the cloud slab 



Name 


Description 


Range 


Distance profile from the clouds to the sensor (km). 
The two clouds are assumed to have the same 
distance. 


Coverage 


The field of view coverage profile of the front cloud 
slab (CI) in the form of a fraction from 0 to 1. 


T c , 


The temperature profile of the front cloud, CI 


NumCompCl 


Number of components in the front cloud 


CompIDi, c i 


The ID tag for the first component in CI 


CLcompl.cl 


The CL profile for the first component in CI 
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The ID tae and CL nrofile for the i th comnonent in 
CI 


CompID num Compl,cl 


The ED tag for the last component in CI 


numComp 1 ,c 1 


The CL nrofile for the last comnonent in CI 




The temnerature nrofile of the hack cloud C2 


NumComp C2 


Number of components in the back cloud 


CompIDi,c2 


The ID tag for the first component in C2 


CL CO mpl,c2 


The CL profile for the first component in C2 




The ID tag and CL profile for the i th component in 
C2 


CompED num Comp 1 ,c2 


The ID tag for the last component in C2 


CL numComp 1 ,c2 


The CL profile for the last component in C2 



The profile as stated in Table 2 may denote the values from 
the starting scan to the ending scan in the simulation run 
(experiment). The profile, which is shown in Figure 10, may be 
captured with 4 numbers: startValue 71, endValue 72, breakPointl 
73, and breakPoint2 74. StartValue 71 does not need to be 
smaller than endValue 72, and the breakpointl 73 may be less 
than or equal to breakPoint2 74. Breakpointsl 73 and 
breakpoint2 74 may be expressed in percentage. Numerous forms 
of profiles, such as constant, linearly increasing, linearly 
decreasing and so forth may be specified with these four 
parameters. Figure 10 shows a profile 75 that may specify the 
values from the start scan to ending scan. 

Another input related to the cloud slabs is the absorption 
coefficients file 15, which may store the absorption 
coefficients of all the agents, simulants and interf erents . The 

19 
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format of absorption coefficient file 15 may be defined as 
having an open-ended data file, fixed- length records; one agent 
per record, and vice versa; and an agent "name" in user 
interface 63 linked to stored index/key. 
5 A given record may contain four elements: Cstring - 

AgentName; DWORD - IndexKey; CString - Comment; and float - 
AgentCoeff icient [2184] . AgentName may be a unique, recognizable 
identifier, 2 -> 16 characters in length. IndexKey (also 
unique) may be transparent to the user, and assigned at the 

10 start of the respective record. Comment may be optional, 
variable (or limited) length. AgentCoeff icient array may 
consist of 2184 floating point values. Upon retrieval, these 
may be converted to mg/m2 by multiplying stored value * 0.0001. 
A starting wavenumber (each 1st data point) may be 388.17558 cm" 1 

15 and a step may be 0.96716738 cm" 1 . For reference, the ending 
wavenumber thus may calculate out to something like 2499.50197 
cm' 1 . 

As to outputs of the simulation software, most of them may 
be stored back into the CADSE database. These outputs may 
20 include simulated interf erograms 64 as scan representations in 
the CADSE database. These scans may be associated with a new 
global file header and subfile information. Another output may 

20 
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be a new global file header. Since the contents of the scans 
(subfiles) may change, it could be beneficial to have a new 
global file header associated with the simulated data. In 
particular, the following information in the global file header 
5 may be changed. A filename may change from d******* to 

S******* 7 where D stands for data and S stands for simulation. 
The rest of the character sequence in the file name may be the 
same. One may note also that the same interf erogram could be 
used multiple times to simulate different 

10 agents/simulants/interf erents , yielding many interf erograms . 
Thus, the association of the synthesized data with a global 
header may take into account this situation. Another output may 
be subfile information originally associated with simulated 
interf erograms 64. An output may also be simulation ground 

15 truth information associated with each of the simulated 

interf erograms 64, i.e., components in the two cloud slabs 
(their spectral ID, temperature, CL, range, FOV coverage) . 
Other outputs that are not stored in the CADSE database may 
include all user inputs to the simulation which need to be 

20 stored in a file with unique name. 

A high level operational flow diagram of the simulation 
software is shown in Figure 11. A user interface and data 

21 
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retrieval block 81 may have an output connected to an input of a 
MODTRAN atmospheric simulation block 82, a sensor self radiation 
and response block 83 and a cloud radiance and transmittance 
block 84, respectively. Another output of user interface and 
5 data retrieval block 81 may go to an input of a synthesis block 

85. An output of MODTRAN atmospheric simulation block 82 may go 
to an input of synthesis block 85. An output of sensor self 
radiation and response block 83 may go to an input of synthesis 
block 85. An output of cloud radiance and transmittance block 

10 84 may go to an input of synthesis block 85. An output of 

synthesis block 85 may go to an input of a data storage block 

86. Sensor self radiation and response block 83, cloud radiance 
and transmittance block 84 and synthesis block 85 may be 
contained in a Matlab module 87. 

15 User interface and data retrieval block 81 v may provide a 

user an interface such that the user can input the parameter 
values. The user interface may be in a form of a sequence of 
windows. Based on the input values, this block 81 may then 
retrieve and compute all of the required data for the 

20 simulation, e.g., the ambient blackbody, and liquid nitrogen 

references. Four categories of outputs defined for this block 
81 may include an input file to MODTRAN block 82, inputs to 

22 
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sensor response/self radiation computation block 83, inputs to 
cloud radiance and transmission computation block 83, and inputs 
to synthesis block 85. 

MODTRAN simulation block 82 may compute the foreground 
atmospheric radiance and transmission. Input to block 82 may be 
in a form of an input file. Outputs of block 82 may also in the 
form of files. Decimation of the MODTRAN outputs to match the 
simulated resolution may be required. 

Sensor response and self radiation block 83 may compute the 
sensor self radiation and the sensor response, which may be used 
for calibration and uncalibration . Calibration may remove the 
sensor characteristics and artifacts before the 
agent/interf erent spectral characteristics are added. 
Uncalibration, then, may add back the sensor characteristics and 
artifacts to the simulated spectrum. Computing the sensor self- 
radiation and the sensor response may require ambient black body 
reference interf erograms , liquid nitrogen (LN 2 ) reference 
interf erograms , and ambient temperature. The interf erogram to 
spectrum transformation may be a commonly used function that 
involves some computations. The outputs of this transformation 
may be the phase corrected spectrum and the phase array. 
Likewise, the inverse transformation of a spectrum to an 
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interf erogram may be another commonly used function. The output 
of this transformation may be the interf erogram. 

Cloud radiance and transmission block 84 may compute the 
transmission and radiance of the two cloud slabs based on their 
5 temperatures, their components' absorption coefficients, and 
CLs . The computations for each of the cloud slabs may include 
computing the exponent of the transmission based on the 
components and their absorption coefficients (ARG as described 
in the simulation model) , computing the transmission from the 

10 exponent, computing the radiance from the transmission and 

theoretical black body at the cloud temperature, and resampling 
the transmissivity and radiance to the synthesized signature's 
sample comb. The spectral range of interest may be from 700 cm" 1 
to 1450 cm" 1 . However, the simulation may need to extend this 

15 range further to cover any filter window, and ensure that 

discontinuity, which causes undesirable disturbance in the final 
signal, is avoided. 

Synthesis block 85 may synthesize the agent /interf erent 
into the input interf erogram. Computations may include coadding 

20 (average) the number of interf erograms based on the numCoadd 

input value, subtracting averaged liquid nitrogen interf erogram 
from the averaged interf erogram, transforming the difference 
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interf erogram to a background spectrum (I B -l) / computing the 
calibrated spectrum, or the scene radiance (N B ) , synthesizing the 
simulated spectrum (N SIMUL ) , computing the sensed simulated 
spectrum (I S imul) , adding appropriate noise to the simulated 
spectrum (I ' simul) / computing the inverse transformation to obtain 
the synthesized interf erogram, and adding averaged liquid 
nitrogen interf erogram in order to negate the effects of 
normalization for the simulated interf erogram. 

Data Storage block 86 may store the results into the CADSE 
database and to the output file. Modifications of appropriate 
attributes and creation of new attributes in the database may be 
necessary. Interface between the Mat lab and the CADSE 
simulation module may be specified in the following. CADSE may 
invoke the Matlab function and pass the corresponding arguments. 
The Matlab function may be called Main_Simulat ion as described 
below: 

function [jslscad_interferograni] = Main_Simulation( subfiles, sf_specs, sf_A, sf_L, 

sf_calib_file_specs, configuration_specs, 
absorptivitycoefficients, N_F, tauJF, step_size, ilsf, 
halfwidth) 

where 

Inputs may include: 

subfiles (NPTS x nscan) matrix for storage of background interferogram 

sf_specs nScan x 3 matrix for storage of sensor gains (col 1), centerburst 

(col 2),calibration file offset(col 3) 
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sf_A NPTS x nCalibration matrix for storage of ambient blackbody 

calibration interferograms 

sfJL NPTS x nCalibration matrix for storage of liquid nitrogen 

calibration interferograms 
5, sf_calib_file_specs nCalibration x 5 matrix for storage of sensor gains for ambient 

blackbody calibration interferograms, centerburst locations for 
ambient BB calibration interferograms, sensor gains for liquid 
nitrogen calibration interferograms and centerburst locations for 
liquid nitrogen calibration interferograms respectively 
10 configuration_specs nscan x 23 matrix for storage of the configuration specifications 

per definition shown in TABLE specs below. 

absorptivity_coefficients nrl x nAgent matrix for storage of agent absorptivity 

coefficients 

N_F nrl x nScan matrix for storage of atmospheric spectral radiance 

15 tau_F nrl x nScan matrix for storage of atmospheric spectral 

transmissivities 
step_size interferogram stepsize [micrometers] 

ilsf instrument line shape function 

half_width half width of time domain weighting function for calculation of 

20 reference spectrum. 

Outputs may include: 

jslscadinterferogram NPTS x nScan array for storage of synthetic interferogram 
25 The number of points per interferogram, NPTS, may be defined internally. 

TABLE specs may include: 
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Matrix Column 


Definition 


1 


Ambient Temperature 


2 


Spectral Resolution 


3 


FOV coverage 


4 


Cloud Slab #1 - temperature 


5 


Cloud Slab #1 - number of components 


6 


Cloud Slab #1 - component 1 absorptivity coefficients 


7 


Cloud Slab #1 - component 2 absorptivity coefficients 


8 


Cloud Slab #1 - component 3 absorptivity coefficients 


9 


Cloud Slab #1 - component 4 absorptivity coefficients 


10 


Cloud Slab #1 - component 1 CL 


11 


Cloud Slab #1 - component 2 CL 


12 


Cloud Slab #1 - component 3 CL 


13 


Cloud Slab #1 - component 4 CL 


14 


Cloud Slab #2 - temperature 


15 


Cloud Slab #2 - number of components 
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16 Cloud Slab #2 - component 1 absorptivity coefficients offset 

17 Cloud Slab #2 - component 2 absorptivity coefficients offset 

1 8 Cloud Slab #2 - component 3 absorptivity coefficients offset 

19 Cloud Slab #2 - component 4 absorptivity coefficients offset 
5 20 Cloud Slab #2 - component 1 CL 

21 Cloud Slab #2 - component 2 CL 

22 Cloud Slab #2 - component 3 CL 

23 Cloud Slab #2 - component 4 CL 



10 A two-slab cloud model may be used in the agent/interf erent 

simulation software. Figure 2 shows the geometry of the two- 
slab cloud model for the agent /simulant /interf erent simulation. 
The simulation may model the sensed radiance as the cumulative 
effects of the background radiance, the radiance and 

15 transmission of clouds 91 and 92 (CI and C2) and that of the 

foreground. Thus, the four components (i.e., background, Cloud 
91 (CI) , cloud 92 (C2) , and foreground) of the model may have 
their own radiance, N B ' , N C i, N C 2 / and N F , respectively. The two 
cloud slabs and the foreground may have their own transmission, 

20 T C i, t C 2/ and x F , respectively. The two cloud slabs may be assumed 
to be at the same distance from sensor 21. Each cloud (91, 92) 
may be composed of a homogeneous mixture of chemical components 
at a specific temperature. The simulation may allow a maximum 
of four components in each cloud slab. Each component may be 

25 represented by a set of absorption coefficients at a specific 

concentration x path length (CL) . In addition, the back cloud, 
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C2 , may be assumed to have a 100 percent field of view (FOV) 
coverage but the front cloud, CI, may have a 0 to 100 percent 
FOV coverage. A shadow factor, s x , which has a value of 0 to 1, 
may be defined as the fraction of cloud 91 (CI) covering cloud 
92 (C2) within the FOV coverage. Thus, Si = 1 may imply that 
both clouds 91 and 92 are within the FOV, and s x = 0 may imply 
that cloud 91 (CI) is outside of the FOV. For convenience, 
another term, s 2 , may be defined as 1 - Si . 

The model may implicitly simulate one cloud slab. Several 
methods may be used to accomplish this. Setting s x = 0 is one 
way; setting the number of components in one of the clouds to be 
zero is another way. This model may assume that sensor 21 
measured the sensed signal ( interf erogram) without the two 
clouds, I B . The sensed interf erogram, I B , may be converted to 
radiance without the cloud, N B . One may assume that N B = x F N B + 
N F . Then the simulated radiance, N S imul, entering the sensor may 
be 

NsiMUL = Sx (N B -N F ) T C iTc2 + S 2 (N B -N F ) T C 2+N f + SiNc2XciT F +S 2 Nc2Xf + S 1 NciT f ■ • ■ m 

This simulated radiance may passes through sensor 21, which may 
detect the radiance and output the simulated sensed radiance, 
Isimul a.s 

1 ' simul = N SIMUL ,*R 
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This equation does not take into account the artificial 
attenuation of the noise originally in the measured background 
spectrum. (The noise compensation may take place in the raw 
spectrum after the simulated radiance spectrum has been 
5 uncorrected for sensor spectral response.) One may note that 

IsiMUL = 1 ' SIMUL + N no i S e . 

For now, the noise spectrum, N no i se , may be assumed to be 0, 
because the noise is not estimated accurately. A more 
definitive method to capture or compute the noise may be needed. 

10 Isimul may be in the spectral domain. Following computations of a 
method 1 may transform Isimul to an interf erogram: 1) Restore 
phase from the reference phase (computed in the forward 
transformation) ; 2) Normalize the phase compensated measurement 
data plus agent spectrum by the phase compensated measurement 

15 data spectrum, I B ; 3) Perform inverse FFT ; 4) Truncate the 

transformed interf erogram to half its length (assume zero-filled 
before) ; 5) Remove original apodization effect by dividing the 
truncated interf erogram with a window shown in Figure 12 (8 is 
the maximum peak location 93 in the figure) ; 6) Combine the gain 

20 compensated measurement data interf erogram with the unapodized 
difference scaled by a factor of 2.0; and 7) Un-normalize with 
original gain . 
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Following may provide all the equations for computing the 
terms in the above equations. Si is the shadow factor, which may 
be the fraction of cloud CI covering cloud C2 within the FOV. 
Thus, s 2 may equal 1 - Si. 

R may be a response of the sensor and is computed as: 

R= ! A ~ L i 

Nbb-N bb 

where 

jy A = the theoretical black body at ambient temperature, 

and 

jsf L BB = the theoretical black body at the temperature of 
liquid nitrogen (77° K) . 

The theoretical black body may be computed using Planck's 
function with 

n * bb = cix o 3 I (e c2 "/ T - i) w/cmVcrrfVst / 

where 

Ci = 1.19089 x 10" 12 w cm 2 sr" 1 , 

c 2 = 1.4388 cm • °K, 

o = wave number value (cm" 1 ) , 
T = temperature (°K) , and 
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N B = total radiance of the background and foreground, 
(without the cloud slabs) . 
N B may be computed from calibrating the spectrum, I b -l, which is 
the measured background signal with the sensor self radiation 
removed, as 

N B = 
B R 

I b -l and I A -l may be computed from the background (B) and the 
averaged ambient blackbody (A) interf erograms (referred to as 
the input interf erogram below) respectively. They may be 
computed with a method 2 as follows: 1) Normalize the input 
interf erogram for gain; 2) Subtract the normalized averaged 
liquid nitrogen interf erogram from the normalized input 
interf erogram; 3) Apodize the difference interf erogram using a 
triangular apodization function, which has a value of 0.5 at the 
center burst location (location of zero phase reference) and 0 
values at the extended points 0 and Npts+1; 4) Zero-fill the 
apodized interf erogram to double the length; 5) Apply FFT to the 
zero-filled interf erogram; and 6) Phase correct the resulting 
spectrum with the generalized version of the Mertz method that 
accounts for both magnitude and phase. 

One may note the following term definitions. 
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N F/ x F = the radiance and transmission of the foreground; 
N C i, t C i = the radiance and transmission of the front cloud 
slab, Cl; and 

N C 2, t C 2 = the radiance and transmission of the back cloud 
slab, C2 . 

The terms may be computed as 
r cl = e ~ ARG cy an d r C2 = e~ ARGc2 , 

N Cl= N C BB^- T Cx) and N C2 = N C BB^^ T Cl) > 

where 

ARG C( = a X CL , + a 2 CL 2 + ... + a n CL „ / 
ARG C2 = a x CL x + a 2 CL 2 + ... + a m CL m , 

OLi = absorption coefficient for component i as a function of 
o , 

and 

TVab = T ^ e theoretical blackbody at the temperature of cloud 
i , for i = 1,2. 

The following function, f (x) , may be used in the decimation 

process -- that is to convolve it with a high resolution sampled 

sequence, followed by an interpolation and subsample. The 

resulting sequence may have a lower resolution. 

The convolution function may be 
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f(x) = sinc 2 Jc = 



sin 2 x 



(for triangular apodization) , and 



x = 0.88 7i v I f whm, 



where 



u = displacement from center of function (cm 1 ) , and 



5 



fwhm = full width at half maximum. 



For JSLSCAD-like data, one may use a 3 3 point function, 
normalized to unit area. If the high resolution sample interval 
is 1.84451968, and the laser diode frequency is 0.662, then 

o = 1.84451968 x INDEX (from center of function) cm" 1 , and 



It may be noted that the above discussion is in the 
spectral domain. The input to the simulation may be an 
interf erogram. Thus, the simulation may first transform the 
input from the time domain to the spectral domain. Then the 

15 agent and interf erent synthesis as outlined may be performed. 
Afterwards, the inverse transformation may convert the 
synthesized spectrum back to the simulated interf erogram. 

Although the invention has been described with respect to 
at least one illustrative embodiment, many variations and 

20 modifications will become apparent to those skilled in the art 
upon reading the present specification. It is therefore the 



10 



fwhm = 3.2. 
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intention that the appended claims be interpreted as broadly as 
possible in view of the prior art to include all such variations 
and modifications . 
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